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Quasars are the visible signatures of super-massive 
black holes in the centres of distant galaxies. It has 
been suggested^ that quasars are formed during "ma- 
jor merger events" when two massive galaxies col- 
lide and merge, leading to the prediction that quasars 
should be found in the centres of the regions of largest 
overdensity in the early Universe. In dark matter 
(DM)-dominated models of the early Universe, mas- 
sive DM halos are predicted to attract the surround- 
ing gas, which falls towards its centre. The neutral gas 
is not detectable in emission by itself, but gas falling 
into the ionizing cone of such a quasar will glow in 
the Lyman-a line of hydrogen, effectively imaging the 
DM halo^. Here we present a Lya image of a DM halo 
at redshift 3, along with a two-dimensional spectrum 
of the gaseous halo. Our observations are best un- 
derstood in the context of the standard model for DM 
halos^; we infer a mass of (2 — 7) x 10^^ solar masses 
(Mo) for the halo. 

Using radiative transfer calculations, Haiman & 
Rees^^^^ predicted that gas falling into a DM halo between 
redshifts 3 and 8 that was harbouring a quasar should 
be detectable in Lya-emission at flux levels accessible 
to present day telescopes, owing to the reprocessing of 
quasar ultraviolet photons. Recently, Barkana & Loeb^ 
found absorption features in quasar spectra, which they 
interpreted as a signature of neutral hydrogen (Hi) falling 
into the DM halos surrounding two quasars. Such absorp- 
tion features can be used to study the gas in one dimen- 
sion (ID) along the line of sight, whereas detection of ex- 
tended Lya-emission can be used in a 2D study, providing 
more constraints on the interplay between gas, quasar ra- 
diation and DM halo. 

In a deep Lya narrow-band image we detected asym- 
metric extended emission North-East of the z = 3 radio- 
quiet quasar Q 1205-30, making this quasar a prime can- 
didate for a study of its DM halo. The extended Lya 
emission was first thought to be related to a foreground 
absorber^, but deep follow-up spectroscopy obtained with 
the FORSl instrument on the ESO Very Large Telescope 
allowed us to measure precise redshifts of quasar, ab- 
sorber, and extended emission, clearly linking the ex- 
tended emission to the quasar, not the absorber. The 2D 
spectrum is presented in Fig. [I] (details of the data reduc- 



tion will be presented elsewhere). In the 2D spectrum the 
wavelength increases along the abscissa, and the position 
on the sky changes along the ordinate. In order to re- 
veal the underlying extended emission we have optimally 
extracted the 2D quasar spectrum by fitting and remov- 
ing the quasar spectral point- spread function^. Spectro- 
scopic detection of kinematically resolved extended Lya- 
emission is not uncommon^' ^, but it does not provide suf- 
ficient information to distinguish between scenarios of - 
for example - nearly edge-on disks and gaseous halos. 
For this, morphological information is needed. In Fig.|2t 
we present a 10 x 10 arcsec^ narrow-band image of the 
extended emission. 

A quasar emits its radiation in an ionizing cone with an 
opening angle of ^, while we observe the whole system 
under an inclination angle of 6 (see Fig.[3l). The surface 
brightness at each point in the vicinity of the quasar is 
calculated by integration of the volume emission along 
the line of sight^. We take the gas density profile to be 
a power-law, nHi(r) = nm,i {r/l kpc)~^, with a slope 
a and a neutral hydrogen density at a distance of 1 kpc 
^Hi,i, because this appears in numerical simulations of 
DM halos to be a good approximation at small radii^^. If 
we assume a mass of the DM halo, we may calculate the 
observed infall velocity in the model. At a given projected 
distance from the quasar the observed velocity is the av- 
erage projected velocity weighted with the emissivity of 
the gas along the line of sight. 

Applying this model to Q 1205-30, our calculation 
shows that the Lya halo should be detectable with 8-m 
class telescopes up to several arcsec from the quasar, as- 
suming a modest amount of neutral hydrogen infall. Con- 
sidering - for now - the case where only one cone is vis- 
ible, the extended emission may be symmetric (if ^ ~ 0, 
that is, the cone is seen close to end-on) or highly asym- 
metric (if ^ ^/2). By fitting the calculated surface 
brightness maps to the observed narrow-band image we 
establish a relation between the opening angle and the 
best-fitting inclination angle. Thus, we are left with the 
opening angle, ^, the slope, a, and the neutral hydrogen 
density at 1 kpc, nHi,i, as the only free parameters of the 
model. In Fig.|2l we compare our observed narrow-band 
image to calculated Lya surface brightness maps with an 
opening angle ^ = 110° and three different inclination 
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angles: the best-fit and two others. For the expected large 
opening angles ^ 90° — 120°)^^'^^, varying the open- 
ing angle primarily affects the shape of the surface bright- 
ness profile seen on the side where the emission is weakest 
(in Fig. 13 that would be on side A of the sightline), and 
only has a minor effect on the shape of the main emission 
profile (on side B). The gas density scale, nHi,i, deter- 
mines the normalization of the surface brightness profile, 
but not the shape, which is set by the slope a. The results 
are not dependent on the choice of a power-law gas den- 
sity profile. An exponential gas density profile provides 
an equally good fit to the data. 

The extended emission is spatially resolved, so we can 
measure its surface brightness profile and velocity profile 
(see Fig. |4|). The surface brightness profile is measured 
by integrating the flux from 4900 A to 4947 A in each 
spatial bin. The velocity in each spatial bin is measured by 
fitting a gaussian to the line profile. The infall velocity is 
calculated relative to the redshift z = 3.041 of the quasar. 

We find that the model describes the observations well. 
The opening angle is the only parameter which remains 
unconstrained, so we plot the surface brightness profile 
for a range of opening angles (Fig.|4t). The best-fit val- 
ues for the neutral hydrogen density at 1 kpc vary less 
than a factor two between ^ = 90° and 170°, while the 
slope varies between 0.02 and 0.09. The velocity profile 
of the extended emission (Fig. |4|3) arises as a projection 
effect, and fitting velocity curves of canonical DM halo 
profiles^®^^ to the observed velocities, we infer a virial 
mass of (2 - 7) X lO^^M© for the DM halo. An identical 
mass estimate is obtained when using an exponential gas 
density profile. The Hi density from our best fitting model 
allows us to calculate an accretion rate of ~ O.lM0yr~^ 
in neutral hydrogen. The total accretion rate will be higher 
as the gas is highly ionized. 

Where studies based only on spectra are unable to pro- 
vide evidence to distinguish between several different sce- 
narios, the combination of deep imaging and deep spec- 
troscopy allows us to rule out alternative explanations of 
the extended emission. Jets are generally believed to be 
present in radio-quiet quasars and are predicted to extend 
out to ~ 0.1 kpc^^, where the extended hya emission 
around Q 1205-30 extends out to ~ 30 kpc. Outflowing 
galactic winds are believed to be triggered by supernovae 
going off inside the galaxy, and are therefore expected to 
be metal-enriched. Our detection limit of 4 x 10~^^ erg 
s~^ cm~^ arcsec"^ (3a) ensures that we are able to de- 
tect the Civ and Hell lines typical in extended emission 
around radio-loud quasars We did not detect any of 
these lines, so the gas appears not to have been enriched 
by supernovae. 
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Figure 1: Two-dimensional spectrum of Q1205-30 and the ex- 
tended Lya emission. The wavelength increases along the ab- 
scissa, the position on the sky changes along the ordinate. The 
spectroscopy was performed with a slit position angle of 7.9° 
east of north (see Fig. |2^). a, The 2D quasar spectrum. The 
extended Lya emission is faintly visible at 4, 920 A. b, The 
quasar spectrum has been subtracted between the wavelengths 
4, 850 — 5, 000 A, clearly revealing the underlying extended 



Figure 2: Narrow-band image of the extended Lya emission 
compared to models, a, Image of the 10 x 10 arcsec^ region 
around Q 1205-30 seen in Lya. The quasar has been subtracted 
from this image to reveal the underlying extended emission. A 
1 .2-arcsec-wide slit is overplotted at a position angle of 7.9° east 
of north, and the position of an unrelated foreground galaxy^ is 
marked gl. The vertical bar on the lower left is 1 arcsec high, 
b - d. Calculated 10 x 10 arcsec^ surface brightness maps for 
an opening angle of ^ = 110° and three inclination angles. 
The best fitting inclination angle (6 — 27°) is shown in b, two 
other inclination angles are shown in c = 10°, excluded at 
3.5a) and d = 50°, excluded at 3a) for comparison. Given 
an opening angle, the best-fit inclination angle is determined in 
the following way. An average surface brightness profile, paral- 
lel and perpendicular to the apparent major axis of the extended 
emission, is fitted to the same average calculated from the mod- 
els. For each opening angle ^ the best-fitting inclination angle 
is found. 
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Figure 3: Schematic presentation of the model, a, A quasar 
is located in the centre of a DM halo into which neutral hydro- 
gen is falling. The ionizing photons from the quasar are emitted 
in a cone of opening angle ^, and they cause the infalling Hi 
to glow in Lya. The observer views the system inclined at an 
angle 0. The surface brightness at each point is calculated by 
integration of the volume emission along the line of sight^^^^: 
SLya = ^^^"^^^^^^^ where Ei^ya. = 10.4 eV is the 
energy of a hya photon, nLya {r) is the production rate of hya 
photons at a distance r from the quasar, is the quasar lumi- 
nosity distance, and Da is the angular distance, such that ^ is 

the conversion from cm^ to arcsec^ . The hya production rate is 
nLya(r) = rithinnm{r)r{r), and ?7thin = 0.42 is the probabil- 
ity for an ionizing photon to result in a Lya photon in the opti- 
cally thin case, nm {r) is the volume density of neutral hydrogen 
atoms, and r(r) is the ionization rate of hydrogen atoms, b, A 
schematic illustration of the predicted surface brightness pro- 
file for various opening angles ^ = 90° (dotted), 110° (solid), 
130° (long-dashed), 150° (solid), and 170° (short-dashed), and 
a fixed inclination angle = 40°. The profiles are presented 
as they would appear when observed with a seeing of 0.7 arc- 
sec through a slit aligned with the line between the sightline / 
and the symmetry axis of the cone. Negative angular distances 
correspond to emission seen from side A, and positive angular 
distances correspond to emission seen from side B. 



Figure 4: Calculated surface brightness profile and projected 
velocity, a. Comparison of the calculated and observed sur- 
face brightness profile. Thin lines are the calculated profile for 
various opening angles, and data points are our observations of 
the extended Lya-emission around Q 1205-30. The calculations 
were performed for opening angles and inclination angles of 
^ = 90°, e = 22.5° (dotted), ^ = 110° , = 27.2° (solid), 
^ = 130°, = 31.8° (long-dashed), ^ = 150°, = 36.5° 
(solid), and ^ = 170°, = 41.1° (short-dashed). The thick 
dashed line shows our 5cr detection limit, b. Comparison of cal- 
culated and observed projected velocity relative to the quasar 
redshift. Lines are the calculated best-fit velocity profiles for 
opening angles and inclination angles of ^ = 90°, = 22.5° 
(dotted), ^ = 110°, = 27.2° (solid), ^ = 130°, = 31.8° 
(long-dashed), ^ = 150°, = 36.5° (solid), and ^ = 170°, 
6 = 41.1° (short-dashed). Three of the measurements in panel 
a were too faint to allow a secure velocity determination. In 
both panels the angular distance was converted to physical dis- 
tance assuming a flat Qa = 0.7 Universe with a Hubble constant 
Ho = 70 km s~^Mpc~^. Error bars are la. 



